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Abstract. We show how two qubits encoded in the orbital states of two quantum 
dots can be entangled or disentangled in a controlled way through their interaction 
with a weak electron current. The transmission/reflection spectrum of each scattered 
electron, acting as an entanglement mediator between the dots, shows a signature of the 
dot-dot entangled state. Strikingly, while few scattered carriers produce decoherence 
of the whole two-dots system, a larger number of electrons injected from one lead with 
proper energy is able to recover its quantum coherence. Our numerical simulations 
are based on a real-space solution of the three-particle Schrodinger equation with open 
boundaries. The computed transmission amplitudes are inserted in the analytical 
expression of the system density matrix in order to evaluate the entanglement. 
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1. Introduction 

Among the various proposals and schemes advanced for realiable quantum computing 
architectures [TJ [2j [3j H] , semiconductor double quantum dots (DQDs) are considered 
very promising candidates for the realization of quantum bits and gates [31 IH El El E] ■ 
Indeed, these structures can arbitrarily be scaled to large systems and could be easily 
integrable with other microlectronic devices. Besides their potentialities in the frame 
of quantum information science, DQDs are also very interesting from the basic physics 
point of view, as they enable both to analyze the peculiar features of electron transport 
phenomena, and to relate them to the appearance of quantum correlations [HI El [101 E] • 
A number of implementations of DQDs qubits have been investigated from the 
theoretical and experimental points of view [31 HI El El U\- Two degrees of freedom, 
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spin and charge, can be used to encode the qubit. While the feasibility of quantum 
logic gates acting on spin states is hampered by the need of local magnetic fields |T2] , 
state-of-the-art nanofabrication technology p3J Q31 [15] allows for a precise control of the 
local charge and orbital degrees of freedom. Recently, Shinkay et al. [15] have realized 
the coherent manipulation of charge states in two spatially separated DQDs integrated 
in a GaAs/AlGaAs heterostructure. Specifically, multiple two-qubit operations, such as 
the controlled-rotation and the swap, have been successfully implemented. 

The main threat to the correct functioning of quantum information processing 
devices is represented by the decoherence stemming from the interaction with the 
external environment or, from a different perspective, by the uncontrolled entanglement 
of the qubits with the environment. For charge states in semiconductor DQDs, the 
loss of coherence is mainly due to the coupling of the carriers to the crystal lattice 
vibrations and to the Coulomb interaction with other charged particles [IB] . In fact, the 
decoherence induced by the electron-phonon interaction has been widely investigated in 
the literature [T71 UHl HH ED], where the quantum dots (QDs) are usually considered 
as two point-like systems with two energy levels coupled to the phonon bath. Analytical 
estimations of the decoherence effects and of their characteristic timescales have been 
given by using various techniques, ranging from the Born-Markov approximation [IT] to 
perturbative calculations in a non-Markovian regime [7J [18] . On the other hand, the role 
played by the Coulomb interaction between charge carriers into the loss of coherence 
has not been deeply analyzed. 

In this paper, we intend to investigate the entanglement properties of bound 
electrons in a GaAs DQD when other electron pass through the structure. In particular, 
we focus on the appearance of quantum correlations in a three-particle scattering, where 
charged particles incoming from a lead enter, one at a time, in a DQD structure and 
interact via Coulomb potential with two electrons, one in each dot. In our scheme, the 
electrons crossing the device have the double role of entanglement "mediator" between 
the two dots and between the DQDs system and the leads, i.e. the environment. 
The aim of this work is to show how the Coulomb interaction between the system 
and the "mediator" can be, under certain conditions, a suitable means to entangle or 
disentangle in a controlled way the qubit states encoded in the single-particle energy 
levels of the dots. A detailed theoretical estimation of such effects results to be of 
interest also for the experimental feasibility of the above qubit, as the results reached in 
the production, manipulation, and coherent control of charge states in DQDs seem to 
indicate [131 US] • I n this view the entanglement/disentanglement of the DQDs can be 
connected to the precise engineering or the suitable tuning of physical and geometrical 
parameters, such as the DQDs level spacing or the current intensity, here modelled as 
the successive injection of carriers in the scattering region. 

The numerical procedure used to solve the model is a generalization of the quantum 
transmitting boundary method (QTBM) [211 [22J. It allows us to find the reflection and 
transmission amplitudes of each scattering channel as a function of the initial kinetic 
energy of the incoming carrier and of the state occupation of the dots. Our analysis is 
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time-independent, in the sense that the few-particle scattering states are obtained by 
the solution of a time-independent open-boundary Schrddinger equation, and therefore 
does not permit to evaluate the dynamics and the characteristic timescales of quantum 
correlations [22j [23] . On the other hand, the QTBM takes explicitely into account the 
spatial structure and therefore the size and the shape of the dots, thus permitting to 
overcome the approximations implied in the description of a dot in terms of a two-level 
point-like system. We single-out the peculiar mechanisms of electron transport through 
DQDs resulting in resonances in the transmission and reflection spectra and thus leading 
to entanglement and decoherence [2U [251 1231 [261 121] • In the evaluation of such effects, 
both the transmitted and reflected components of the scattered wavefunction are taken 
into account. 

The paper is organized as follows. In Section [21 we introduce the physical system 
reproducing a DQD structure in GaAs, and illustrate the computational approach 
adopted to find the few-particle scattering states. The description of the DQD in 
terms of a two-qubit model and the discussion about the theoretical procedures used 
to evaluate the entanglement and the decoherence are given in Section [3j In Section [H 
we show the numerical results first obtained for the scattering of a single carrier, then 
for a weak electric current. For the latter case, the condition leading to the maximum 
entanglement production or to the complete disentanglement are analyzed in detail. 
Finally, in Section [5] we comment on the results and draw our conclusions. 

2. The DQD structure and the computational approach 

In our model (see Figure [1]), we consider an electron incoming from the left, with respect 
to the DQD one dimensional (ID) structure, with kinetic energy To. We examine the 
case of one scattered carrier at a time, i.e. we suppose that an electron enters the 
scattering region only after the previous one has already left. Such an assumption, and 
the fact that the charging energy of the DQD is larger than the spacing between the two- 
particle energy levels, means that our system always operates in three-particle regime. 
The incident particle is scattered via Coulomb interaction by the two electrons bound in 
a structure potential V s in a region of length L (Figure [1]) mimicking a DQD structure 
and constituted by two potential wells separated by a potential barrier wide enough to 
make negligible the Coulomb interaction between the two confined particles. Moreover, 
the structure is connected to external leads kept at zero potential. The N two-particle 
bound states and energies of the DQD will be indicated by |S n ) and e n , respectively 
(with nO indicating the ground state). As it will be shown in the following, some of them 
can also be expressed in terms of \xfXm) an d (Ep 1 + E^), where |xf )(\Xm)) indicates 
the single-particle bound state of the right (left) dot with energy Ep(E^). Due to the 
symmetry of the potential V s , Ej 1 = Ep. 

As anticipated, we restrict our investigation to a ID analysis. Such an assumption, 
needed to solve numerically the few-particle problem with open boundaries, is physically 
reasonable if the transverse dimension of the structure is small compared with the other 
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Figure 1. Profile of the potential V s (x) in the scattering region of length L=100 nm: 
the two potential wells are 110 meV deep and 30 nm wide and are separated by a 20 nm 
barrier. The dashed lines indicate the single-particle energy levels Eq, Ei, and E2 of 
the ground, the first and second excited state of the dots respectively, with _Eo=-105.6 
meV,£a=-92.5, and i?2=-71.5 meV. With our parameters the levels E2 (and above) are 
found to have negligible occupancy. In our numerical calculations we take m*=0.067 
m e , with m e indicating the bare mass electron, and £=12.9. 



length scales. In this case, all the particles can be supposed to occupy the lowest 
single-particle transverse subband [28] ■ Furthermore, we consider for the incoming 
electrons only energies below the ionization threshold of the DQD. This means that 
when the outgoing electrons leave the scattering region, either reflected or transmitted, 
the confined particles remain in a bound state of the DQD. 
The three-particle Hamiltonian H is given by 

H(xi,x 2 ,x 3 ) = , Ho(x 1 )+'Ho(x2)+'Ho(x 3 )+^2^2 exp (-ry/A d ), (1) 

t=i j=i 



where = y (x{ — Xj) + d 2 and Tio(xi) is the single-particle Hamiltonian 

e and m* are the dielectric constant and effective mass of the GaAs, respectively. The 
term describing the mutual interaction between particles in the r.h.s. of equation ([1]) 
is a screened Coulomb potential with a Debye length A^, here taken significantly larger 
than the characteristics length of the structure. Furthermore, the Coulomb term also 
accounts for the transversal dimension d of the confined system (with d—1 nm) through 
a cut-off term. In our approach, the fermionic nature of the carriers is explicitly 
accounted for by antisymmetrizing the two- and three-particle wavefunctions S n (x 2 ,x 3 ) 
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and ty(xx,X2, x 3 ) for any two-particle exchange. It is worth noting that the Hamiltonian 
given in equation ([1]) does not include spin-orbit terms. As a consequence, since the 
orbital wavefunction is antisymmetric, we are simulating a three-particle system with 
a symmetric spin component, as in the case of three spin- up (or spin-down) electrons. 
Furthermore, we do not include electron-phonon interaction. Infact, the aim of our work 
is to investigate the role played by the Coulomb interaction among the electrons in the 
system into the appearance of entanglement and decoherence. Our device is supposed 
to operate at a temperature in mK regime, where only spontaneous emission is effective 
and ignoring the coupling of the electrons with the surrounding crystal lattice does not 
constitute a crucial approximation, as we will describe in the following. In GaAs-based 
structures with levels splitting of few meV, single-phonon processes lead to an excited- 
state lifetime of the order of 10 -10 s, while multi-phonon and multi-electron processes 
are orders of magnitude less frequent (see References [29j [30]). Given a kinetic energy 
of 15 meV, a single electron traverse our 100 nm long device in about 10 -13 s. If we 
suppose independent electrons injected from the lead at a mean rate of one every 10~ 12 
s, corresponding to a current of about 0.16 fiA, 100 carriers can be scattered through the 
double-dot before a phonon-induced relaxation takes place. As a consequence, in the 
following sections where we consider a weak electric current, we will limit our calculation 
up to 60 electrons. 

The scattering states of the three particles are obtained by solving the time- 
independent open-boundary Schrodinger equation = E^> in the cubic domain 
{xi,X2,X3} with Xi G [0,L] and with the Hamiltonian H given in equation (CQ). For this 
purpose, we applied a few-particle generalization of the so-called QTBM [21], allowing 
one to include proper open-boundary conditions for each edge of the domain. These 
describe the particle coming from the left as a plane wave with energy T and wavevector 
ko while the other two electrons are set in a two-particle bound state of the DQD 
with energy ej. Moreover, to account for the exchange symmetry of the three-particle 
wavefunction, also the antisymmetry of the boundary conditions is imposed, as shown 
in previous works [221 E3] . 

The correlated scattering state when the particle 1 is localized in the left lead (that 
is x\ < 0) reads 

^(x 1 ,x 2 ,x 3 ) = 

(xi<0) 

M 00 

S^xa, x 3 )e lh0X1 + b jn Z n (x 2 , x 3 )e- ik ^-^ + ^ b jn E n (x 2 , x 3 )e k ^^ (3) 

n=0 n=M+l 

where j = ... iV is the index of the initial two-particle DQD state. Analogously, 
when the particles 2 and 3 are in the left lead, the boundary conditions are 
*(xi,x 2 ,x 3 )|( a:2<0 )=-*(x2,xi,X3)|( a:2< o) and ^(xi, x 2 , x 3 )\ {x3<0) = ty(x 3 ,x 1 ,x 2 )\( X3< y, 
respectively. In the above expression k^_ n )=^2m*T^_ n ~ ) , where T^_ n ) denotes the 
kinetic energy of an electron freely propagating in the lead, as obtained by energy 
conservation Tjj- n ) — T + €j — e n . For sake of simplicity, we set h=l here and in the 
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following. 

The first term appearing in the r.h.s of equation ([3]) describes electron 1 incoming 
from the left lead as a plane wave with energy To, while the other electrons are in the 
two-particle bound state |Bj). The second term represents the linear combination of all 
the energetically-allowed possibilities when particle 1 is reflected back as a plane wave 
with wave vector ktj- n ) and the DQD is in the state |H n ). M indicates the number 
of states for which T^_ n ^ > 0. The last term accounts for those states with Ttj- n ) 
negative, which describe particle 1 as an evanescent wave in the left lead. Therefore, 
the coefficients bj n are the transition amplitudes between the initial state Ej(x 2 , x 3 )e lk ° Xl 
and the final state E n (x 2 , x 3 )e lk ^- n ^ Xl when the incoming carrier is reflected. 

If particle 1 is in the right lead, the three-particle wavefunction takes a form similar 
to expression flH]), which describes the outgoing travelling and evanescent modes of the 
electron in the right lead (x\ > L): 

M oo 

V(x 1 ,x 2 ,x 3 ) =Vc, Il H 1! (x 2 ,x 3 )e^-»> I1 + V c jn E n (x 2 ,x 3 )e- k o-^ x \ (4) 

(xi>L) * — ' z — ' 
v ' n=0 n=M+l 

while ^ (x!, x 2 , x 3 )\( X2>L )=-^ (x 2 , Xi, x 3 )\ {x2>L ) for x 2 > L and ty(xi,x 2 , x 3 )\( X3>L )= 
"$?(x 3 ,xi, x 2 )\( X3> l) for x 3 > L. The coefficients Cj n describe the transition amplitudes in 
the n-th channel, i.e. when the bound particles are in the n-th state of the DQD. The 
boundary conditions are given by equations ([3]) and Q with xi=0 and x\=L. They are 
coupled to the Schrodinger equation and discretized by a finite-difference method. A 
system of seven equations is obtained, whose numerical solution provides the unknown 
coefficients bj n and Cj n , and the three-particle wavefunction ^/(xi,x 2 ,x 3 ) in the internal 
points. 



3. Decoherence and entanglement of the two-qubit model 

As a consequence of the scattering, quantum correlations between the single-particle 
energy levels of the bound electrons and the energies of the scattered electron appear. 
They are responsible both for the loss of quantum coherence of the two-particle state 
of the DQD and for the building up of quantum entanglement between the two dots. 
First, we show that under some approximations the DQD system can be reduced to a 
two-qubit model. Then, we describe in detail the theoretical approach used to evaluate 
entanglement and decoherence in such a model. 

Although the fermionic nature of the carriers has been explicitly taken into account, 
as shown in Section [21 in order to solve numerically the physical system, we do not use 
entanglement criteria for identical particles [3U [32} [33] . In fact, the scattered carrier, 
either transmitted or reflected, can be assumed to be far from the scattering region, 
while the bound particles are essentially trapped into two deep potential wells far away 
from each other. So the spatial overlap between the particles results to be negligible. 
Therefore, the position variables can be used to distinguish the particles, while quantum 
correlations are evaluated between scattering channels [3HE3]. By moving from spatial 
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to energy representation for the quantum states, and taking as input state \^j^) = \T €j), 
which describes the carrier incoming from the left lead with kinetic energy T and the 
particles bound with energy €j, the output \§out) reads 



M M 



\®OUT) = ~ h ^\ T 5-nfn) + Cjn\T^ n) e n ) (5) 
n=0 n=0 

where the coefficients bj n are given by bj n = bj n /{^}2!^ =Q {\bj n \ 2 + |c,,- n | 2 ), and the analogous 
expression holds for Cj n , while \T^^\e n ) indicates the state with the carrier reflected 
(transmitted) as a plane wave (with kinetic energy T( n _j) = fc? .^/2m) and the other two 
electrons bound in |H n ) with energy e n . It is worth noting that in the above expression 
we have omitted the reflected and transmitted outgoing evanescent modes, since their 
contribution to the total current is zero and cannot be responsible of any entanglement. 

So far, we considered the case of the injection of a single carrier in the scattering 
region when the two particles trapped in the DQD structure are described by a pure 
state. Let us now examine the injection of a second electron also with kinetic energy 
To in the scattering region, occurring after the exit of the previous one from the DQD 
structure via transmission or reflection. We indicate as a and b the first and the second 
injected electron, respectively. When b enters the DQD, the bound electrons are not in 
a two-particle pure state, since they are coupled to the energy levels of the carrier a, as 
shown by equation (jSJ). Therefore, the scattering between the electron b and the other 
two bound in the DQD, will give the four-particle state 

M M M M 

\®out) = bjrj)nm\T'^_ m }T^_ n ^E m ) + ^ bj n c nrn \T^_ m ^T^_ n ^e m ) 

71=0 m=0 n=0 m=0 

M M M M 

+ Ci^ nm \T{n-m)T^j- n ) e rn) + Cj n C nm \T^ n __ m ^ T^_ n ^e m ) , (6) 

n=0 m=0 n=0 m=0 

where TL^^T^ + e n — e m and T( has the same meaning of TJ but refers to the second 
scattered particle, i.e. electron b. As more electrons are scattered, the output state 
describing the system involves more and more terms corresponding to the various 
scattering channels. For the sake of simplicity, here we only describe the theoretical 
procedure used to calculate entanglement and decoherence in the case of the injection 
of a single carrier. In the subsection 14. 21 the evaluation of decoherence and entanglement 
due to the interaction with a large number of injected carriers (mimicking an electric 
current) will be given for a specific case. 

In order to compute the non-separability degree of the DQD system into the product 
of single-particle states of the dots, that is the dot-dot entanglement, we have to move 
from the description of the the DQD in terms of two-particle bound states to the 
description in terms of single-particle bound states of the left and right dots. Thanks 
to the negligible Coulomb interaction between the two dots of Figure [TJ the four lowest 
two-particle orbital states |e n ) can be written in terms of single-particle states \E^) and 
\E^), of the right and left dot, respectively (Table [TJ. In particular, |e ) is the product 
of the ground states of the two dots, while |ei) and fa) are degenerate, being each a 
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Table 1. The table displays the scalar product (E^E^\ei) for some values of n, m, I. 
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linear superposition of \EqE h ) and \E^Eq), which represent one electron in the ground 
state of one dot and the other in the first excited state of the other. 1 63) corresponds 
to the first excited states of the two dots. States with two electrons in the same dot 
are included in our calculations but have always a negligible occupancy. This means 
that the coefficients bj n , and Cj n with j > 3 or n > 3 on the right hand side of the 
equation (jHJ) are 0. For our numerical calculations, we used the physical parameters of 
the GaAs material and the DQD potential reported in the caption of Figure [TJ 

In fact, we restrict our analysis to energies of the incoming particle which enable 
up to four scattering channels, i.e. the maximum value of M in the expression ((Sj) is 
3. Under these assumptions a system of two qubits is obtained, in the sense that each 
confined particle can be described in terms of two states: the ground Eq R ^ and first 
excited E^ R ^ energy level of the left (right) QD, encoding the \0l(R)) and \1l(r)) states, 
respectively. Thus, the three-particle quantum state of equation (J5J) can be written as: 



\$ovt) = 



b n + b 



'32 



s/2 



3- 



iOlIr) 



b j2 



s/2 



T< x l L Q R ) + b^\T<_ 3 l L l R ) 



+d j0 \T>0 L R ) - 



V2 



C 7 1 — Cj2 



s/2 



(7) 



where \T>_^) = \T>_ 2 ) and \T^_^j- 



\ 1 j-2l 



2 ), deriving from Eq + E R =E± + E R and energy 
conservation has been taken into account. 

The decoherence undergone by the electrons confined in the dots can be interpreted 
in terms of the lack of knowledge of their quantum state due to the interaction with 
the environment, namely the injected carrier [35J. In other words, due to the coupling 
between the energy states stemming from the scattering event, the two-particle DQD 
cannot be described by a pure state anymore but becomes a statistical mixture. A 
good measure of the degree of uncertainty about such a system and therefore of its loss 
of coherence, is given by the von Neumann entropy of the two particle reduced density 
matrix p r , obtained by tracing the three-particle density matrix p 



matrix representation of p r 
reads 



^out)(^out\ over 
: sc 

in the standard basis B = { (0^0^) , 10^1^), \1l0r), 



the degrees of freedom of the scattered carrier [35]. After the first scattering, the 



\1l1r)} 



Pr 













l/3 + | 2 + l7 + | 2 
+7^7- 






(3 + f31+ 1+1 l 

|/3-| 2 + l7-| 2 

M : 



\ 



(8) 



On demand entanglement in DQDs 



9 



where 



a\ = \b j0 \ + \cj \ 
b n ± b j2 



(3± 
7± 



V2 

Cji ± c j2 



\tu\ 



V2 

IM 2 + |5 i3 | 2 . (9) 

In order to obtain the expression (151) , the orthogonality relations between the states of 
the scattered carrier, (T^^T"- ^) = 5ij and (T<\T>) = Vzj, have been used. 
The decoherence £ can be evaluated by means of the von Neumann entropy as: 

£ = — Tr [p r In p r ] = — | a | 2 In | a \ 2 — rj + lni] + — r]_ In r/_ — \co\ 2 In |w| 2 , (10) 



where 



^ = ^(i/3 + r+i7 + i 2 +i/3-r+i7-r 

±V(\P + \ 2 + I7+I 2 + I/3-I 2 + |7-| 2 ) 2 - 1/3+7- - /3-7+I 2 )- 



It ranges from to In 3/2. For £=0 the two bound particles can be found in a single 
energy level and this implies that no correlation is build up between them and the 
scattered carrier. When the decoherence reaches its maximum, the DQD is found in a 
statistical mixture of the three allowed energies e , e\, and e 3 with equal weight (e 2 is 
equal to ei). This implies that the uncertainty about the system is maximum. 

The reduced density matrix p r describing the bound particles is also used to 
evaluate the dot-dot entanglement through Wootters concurrence C. |36j The latter 
is adopted to quantify the quantum correlations appearing between two qubits which 
cannot be described by a pure two-qubit state because of their coupling with an external 
environment, like in our scenario. C is obtained from the density matrix p r of the two- 
qubit system as [361 EZ] : 



C = max |0, V Ai - V ^2 - V A 3 - V A 4 | , (12) 
where A, are the eigenvalues of the matrix ( = p r (cr^ <g> a^) p* [ay <g> a^) arranged in 
decreasing order. Here (Jy^ is the Pauli matrix ^ Z ^ in the basis {\0l(r)}> \ ^l(r)}}, 

and p* describes the complex conjugation of p r in the standard basis B. The concurrence 
varies from C=0 for a disentangled state to C—l for a maximally entangled state. 

The reduced density matrix p r given in equation (jHJ) shows an X structure, that is 
it contains non-zero elements only along the main diagonal and anti-diagonal. As shown 
in Ref. [38j, for such a class of density matrices the concurrence can be easily evaluated 
and in the case of p r it becomes: 

C = 2max{0, k}, (13) 
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where 



k = |/3+/3* + 7+7* | — = 
(bji + ~ ~ b h) + (cji + c i2 )(c* 1 - ~c) 2 ) 



'jo 



2 + \c 2 



JO 



J j3 



2 +\c 2 



J'3 



• (14) 



C is equal to for |a||o;| > \/3 + ftZ. + 7+7* |, while it reaches its maximum value 1 if 
and only if |/3+/3* + 7+7__|=l and |a| = |w|=0 (and therefore the coefficients bj , c j0 , bj 3 
and Cj'3 vanish). In the latter case, the two-qubit system reduces to a Bell-like state 

l/V^IO^ + e^lUOi^ 



4. Results and discussion 



Here we analyze our numerical results on the decoherence £ undergone by the DQD 
system and the dot-dot entanglement C. We stress again that the former corresponds 
to the entanglement of the DQD electrons with the transmitted/reflected one, while the 
latter is the concurrence between the two bound electrons. In order to single-out the 
specific mechanisms leading to the appearance of quantum correlations, the transmission 
and reflection spectra have been examined. 



4-1- Scattering by a single carrier 

The system has been solved for the potential profile V s (x) sketched in Figure [1] for 
various input states with different energy of the incoming carrier To. In the top- left 
and bottom-left panels of Figure [2] we report the modulus of the transmission (TCs) 
and reflection coefficients (RCs) of the various scattering channels (that correspond the 
energy levels of the two-particle DQDs system) for the input states |T eo) and |T 62), 
respectively, when the incoming electron is injected with a kinetic energy ranging from 
13 to 19 meV. We recall that the two-dot excited state | €2) can be written in terms of 
the qubit states as |e 2 ) = I/v^^JOlIj?) — |1l0_r)^, i.e. it is a Bell state. The TCs and 

RCs are related to the coefficients bj n and dj n given in the expression fl5]). For both 
the input configurations, sharp peaks are present only in small energy interval of the 
reflection and transmission spectra. It is worth noting that for different input states, 
such peaks appear at different values of T . 

To get a better insight into resonances, a zoom of the modulus of TC and RC for the 
various channels in the energy interval around the resonant condition T =15.8 meV is 
displayed in the right panels of the FigureEl While for the input state |Toe2) no resonance 
is present in this range and the sum of the moduli of TC and RC of the channel is equal 
to 1, for |T eo) the spectra show a number of sharp resonances (see the top-right panel of 
Figure [2]). In particular the probabilities of finding the DQD in the state 1 62) (with the 
scattered carrier, either reflected or transmitted) shows a symmetric Lorentzian peak 
around To, while the TCs and RCs of the scattering channel |eo) exhibit a minimum. 
Specifically, the RC presents a symmetric line-shape, and the TC an asymmetric one. 
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Figure 2. Left panels: Modulus of the TC and RC of the channels corresponding to 
the DQD levels |eo) and |ea) as a function of the initial kinetic energy of the incident 
electron T , ranging from 13.2 to 19.2 meV, for two different input states, namely |T eo) 
(top), and |T e 2 ) (bottom): RC of the channel |eo) (solid line), TC of the channel |eo) 
(dashed line), RC of the channel \t2) (dotted line), TC of the channel \e%) (dot-dashed 
line). The modulus of TC and RC corresponding to the other two dot states |ei) and 
| £3} , are zero in the region around the resonance energy Tq= 15.8 meV: for the sake 
of clarity, they have not been reported. Note that the sum of the moduli of the TC 
and RC is 1. Right panels: Modulus of the TC and RC of the channels |eo) and | £2) 
against To close to the resonant energy Tq for the input states |7beo) (top), and |To62) 
(bottom). 



This behavior clearly indicates that in the transmission and reflection spectra different 
kind of resonances appear, namely Breit-Wigner [39J and Fano [40J, respectively which 
have a strong connection into into the building up of quantum correlations, as noted 
elsewhere (231 EJ H2]. The first ones, exhibiting symmetric Lorentzian peaks, stem 
from the coupling of a quasi bound state to the scattering states of the leads, while the 
second ones, characterized by asymmetric lineshapes, are present when two competing 
scattering mechanisms, a resonant one and a nonresonant one, interfere, and are due to 
electron-electron correlation [22| [23]. 

Scattering resonances are a signature of peculiar decohering and entangling 
effects [23], as shown in Figure |3j where the dependence of DQD decoherence £ and 
dot-dot entanglement C are reported as a function of the initial energy of the incoming 
electron around T f° r the input states |T e ) and |T e 2 ). For |T e 2 ), the decoherence 
is practically zero while the entanglement remains 1, as in the initial state. In fact, as 
can be gathered by the behavior of TC and RC, the only coefficients not vanishing in 
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Figure 3. Dot-dot entanglement C and DQD decoherence £ against the initial kinetic 
energy To of the incident electron around the resonance condition To , for two different 
input states: C(|T e )) (solid line), £(|T e » (dashed line), C(|T e 2 )) (dash-dotted 
line), and £(|T e2)) (dotted line). 



equation flH]) are 622 and C22, and the output state \^out) of expression ((7j) reduces to 

\*out) = - (bv\Tf) + £ 2 2|T>)) (\0 L 1 R ) + \1 L R )) . (15) 

This is the factorizable product of a single-particle state of the scattered carrier and a 
two-qubit state describing the bound particles in the DQD, being the latter a Bell state. 
This means that the entanglement between electrons confined in the dots maintains 
its maximum value and also that the quantum information about the DQD system is 
maximal since it is in a pure two-particle state. Thus, the scattering event "preserves" 
the entanglement between the dots while the decoherence effects result to be negligible. 
Such a behavior, evidencing the existence of decoherence-free entangled states of two 
qubits, has been widely discussed in a number of works, which stressed the key role 
played by the symmetric coupling of the qubits with the environment into preserving 
their coherence [I3l EHf H5 l 136] 

When the input state is |T e ), both decoherence £ and entanglement C show a 
maximum where RC and TC of the channel \ez) are resonant (see the top- right panel of 
Figure [2]). In particular £ reaches In 2. Such a value is obtained when the DQD states 
are maximally coupled only to the energy levels T and T_ 2 of the scattered carrier and 
this occurs when the probabilities that the scattering leaves the bound particles in their 
ground |e ) or excited |e 2 ) state are equal. Thus the output state can be written as 

\^out) = (b m \T<) + cool??)) \0 L R )-(b 02 \T< 2 ) + co 2 |T> )) -L (\Q L 1 R ) + \l L Q R )) , (16) 
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where |&oo| 2 + I Coo 1 2 = I ^02 1 2 + | c 2 1 2 ~ 1/2. in this case, from the expressions (TT3~j) and 
( I14p . we observe that the value of the peak of the dot-dot entanglement is equal to 
1/2, as shown in Figure [3j Here an important point to be stressed is that quantum 
correlations are created between the bound electrons, even if their Coulomb interaction 
is negligible due to the large distance between the dots. In fact, even if these can be 
thought as totally decoupled subsystems, the external environment, i.e. the scattered 
carrier, represents the interaction "mediator" and it represents the means to entangle 
them. In the literature, and on the basis of different physical mechanisms, the idea 
of an entanglement mediator has already been used in a number of theoretical and 
experimental models to produce bipartite entangled states [26], H7] SHI HHJ 150] . 



4-2. Scattering by an electric current 

The above results indicate that, for two electrons each bound in the ground state of 
one of the dots, the interaction with a single incident carrier having a suitable kinetic 
energy T , excites the dots. Specifically, the scattering channel corresponding to |e 2 ), 
namely the Bell state describing the first DQD excited level, is activated and quantum 
correlations between the two dots appear, even if the bipartite entanglement production 
does not result to be maximal and immune to decohering effects. Indeed, the probability 
to excite the DQD is smaller than 1. This implies that the two dots cannot be described 
in terms of the state 1 62) alone. Rather, they are in a statistical mixture of ground and 
first-excited states. 

On the other hand, the scattering between a carrier having kinetic energy T and 
the two electrons in the excited maximally-entangled state 1 62) , leaves unchanged the 
DQD state, i.e. the entanglement is preserved and no decoherence effect appears. 
This behavior suggests that the maximum production of entanglement between dots 
set initially in their ground state, can be obtained as a consequence of the successive 
scatterings, one at a time, with carriers injected with energy around T . In fact, at 
each scattering event the probability of finding the DQD system in the excited state | €2) 
becomes larger, and the amount of quantum correlations between the dots increases. 
Such a sequence of carrier injections corresponds to an electric current where all the 
electrons entering in the device have the same energy To. From an experimental point 
of view, such a current can be produced, for example, by using single-electron sources 
such as, electron pumps [9], resonant tunnelling diodes [51], or systems consisting of a 
quantum dot connected to a conductor via a tunnel barrier [52]. All of these mechanisms, 
enable to emit uncorrelated electrons in a given quantum state with a specific energy. 

In order to give a quantitative evaluation of the effect of an electric current on the 

(n) 



DQD state, in Appendix A we have explicitely calculated the reduced density matrix p 



describing the two dots after the injection of n carriers. Its expression in the standard 
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where Poo=|^oo | 2 + |coo| 2 > ranging from to 1, is the probability that a scattering event 
leaves the QDs in the ground energy state when a carrier is injected with kinetic energy 
T . As stated before, for T = T , p 00 is around 1/2. exhibits again an X structure and 
decoherence and entanglement of the system can be evaluated from the equations ( fTOj ) 
and ( 1T4"1) by setting |a| 2 =pQ , /3 + =/3_ = a/ (1 — Poo)/2 an d 7+=7_=w=0. They reads 

e = -^ ln^ -(l-^ )ln(l-p™ ) and C = 1 - p n 00 , (18) 

respectively. When n=0, i.e no scattering occurs, the expression ( !T7l) reduces to 
Pr '' = |0l0r)(0l0r|, which describes the input state where the DQD is in |e ) = 
\0lQr}- For n=l, p-r 1 '' is the reduced density matrix obtained from equation f|T6|) by 
tracing over J$ . In the limit of large n, pr can be written as lim^oo pf 1 = 
| (|0 L lfl)(l fl L | + \0 l 1r)(0r1l\ + \1l0r)(1r0l\ + \1l0r)(0rL l \) which corresponds to 
the Bell state of the two dots ^(\0l1r) + \1l0r)) completely decoupled from the 
environment, with £=0 and C—l. That is, a current of independent electrons (with 
energy T ) entangles the two dots and does not create decoherence. 

Figure H] displays the dependence of the entanglement upon the number n of carriers 
entering in the device at different values of T around the resonant energy T . As shown 
in the inset of Figure HI we find that a series of scatterings does not induce decoherence 
of the DQD first-excited state even for carriers injected with kinetic energies not exactly 
equal but close enough to T . This implies that maximally-entangled states of the DQD 
are produced as an effect of the flux of the charge carriers even if the energy of the 
incident electrons is not precisely the resonant one. Specifically, the farther T is from 
To, the larger n is needed to produce a Bell state decoupled from the environment. 
In fact, when the initial kinetic energy of the carriers gets away from the resonant 
one, the parameter poo, acting as a convergence factor, increases and the injection of 
more carriers in the device is needed to build up the maximum amount of quantum 
correlations between the dots. From the inset of Figure HI we also note that the number 
n of electrons needed to have a vanishing decoherence (i.e. the DQD in a pure state) 
is lower for T closer to T . In particular, £ shows a maximum, whose value is about 
In 2 when the interaction with the injected carriers reduces the state of the system to 
a statistical mixture with equal weights of the ground |e ) and the excited |e 2 ) states. 
This occurs for p l =2~ 1 / n . As expected, the peak is at higher values of n when the 
injection energy of the carriers is farther from Tq and, as a consequence, poo is larger. 

In analogy to the case of entanglement creation described above, a current of charge 
carriers injected at an appropriate energy can be a means to disentangle the DQD 
prepared in the Bell state ^2)^75 (|0l1r) + |1l0#)). In order to show this, Figure [5] 
displays the disentanglement effect for electrons injected with kinetic energy T around 
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Figure 4. Dot-dot entanglement C as a function of the number n of the carriers 
injected into the device evaluated at four different values of the kinetic energy T close 
to the resonance condition To. As input condition, the bound particles occupy the 
ground state of the DQDs system. All the curves tend to 1 for large values of n. The 
farther is the kinetic energy of each incident electron from T , the slower the asymptotic 
value 1 is reached. The inset displays the dependence of the DQD dccoherence £ upon 
n, for the same four values of Tq. 



2.6 meV and scattered by the bound particles of the DQDs. For such a low kinetic energy, 
the scattering by a single carrier leaves unaltered the DQDs system when the bound 
electrons are in the ground state |e ). In fact T is smaller than the energy necessary 
to excite the dots. This means that the the scattered carrier has not been coupled, via 
Coulomb interaction, to the bound particles, which remain maximally disentangled (see 
the bottom inset of the Figure [5]). On the other hand, when the input state of the total 
system is |T e2) the dots can relax. In fact the scattering channels corresponding to |eo) 
show a peak in the transmission and reflection spectra (as shown in the top inset of 
Figure [5j) thus leading to the appearance of decoherence and entanglement. 

By applying the approach adopted above for building up maximum entanglement 
between the bound electrons, we find that the scattering by a current of charge carriers 
with energy around T = 2.6 meV is able to disentangle completely the quantum dots 
without introducing any decoherence. Specifically, after a larger number n of scattered 
carriers, the bound electrons practically occupy the ground state of the DQD system: 
this means C = and £ = 0, as reported in Figure [5j 

5. Conclusions 

The coherent manipulation of electron states is the key ingredient for implementing 
qubits using charge or orbital states of DQD nanostructures. Indeed, it implies the 



On demand entanglement in DQDs 



16 




Figure 5. Dot-dot entanglement C and DQD decoherence £ as a function of the 
number n of the carriers injected into the device at the resonant kinetic energy Tq=2.6 
meV. As input condition, the bound particles occupy the first-excited state |ea) of 
the DQDs system corresponding to the Bell-state -j= (|0j,1.r) + \1l0r))- Both C and 
£ vanish at large n. The top inset displays the modulus of the TC and RC of the 
channels |eo) and | £2) of the DQD as a function of To around the resonant condition 
when the input state of the total system is |To£2): RC of channel |eo) (solid line), 
TC of channel |e ) (dashed line), RC of channel \e 2 ) (dotted line), TC of channel |e 2 ) 
(dot-dashed line). The bottom inset shows the entanglement C and the decoherence 
£ for the two input states |T e ) and |T e 2 )): C(|T e )) (solid line), £(|T e )) (dashed 
line), C(\Tq£2)) (dash-dotted line), and £(|Tbe2)) (dotted line). Note that the abscissa 
scale is the same in the top and bottom insets. 



controlled production or destruction, the manipulation and detection of entanglement 
between the above states. In this spirit, various proposals to produce bipartite entangled 
state have been advanced on the basis of the physical mechanisms requiring two-particle 
scattering, such as the direct interaction between two electrons (2U I25J [231 E31 153] . 
In this work, we have investigated the appearance of quantum correlations between 
the two electrons of a GaAs DQD, as a consequence of the Coulomb scattering by 
one or more charge carriers injected from a lead. We examined the scattering event 
in a three-particle regime (the two electrons trapped in the DQDs and the passing 
carrier, explicitely considered indistinguishable), that is, a carrier is supposed to enter 
in the scattering region only after the previous one has already left. Furthermore, the 
two dots are taken distant enough so that the Coulomb repulsion between the two 
bound electrons is practically negligible. Therefore, unlike other approaches [251 123] . 
the scattered carrier represents the entanglement "mediator", that is, it provides the 
indirect interaction between the particles needed to entangle them. From this point of 
view, various scheme where entanglement between distant particles is produced though 
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their scattering by mobile mediators are present in literature [501 EI]- Unlike our model, 
there the quantum correlations are built among the spin degrees of the freedom of the 
particles. 

A proper tuning of the carrier energy and the DQD geometry reduces the system 
here examined to a simple two-qubit model coupled to the external degrees of freedom 
by the incident electrons. Here, the dots have not been considered as point-like 
systems [7J [T7J HH [191 120] but their effective spatial dimensions are explicitely taken into 
account in the calculation. Indeed, the knowledge of the electron spatial wavefunctions 
corresponding to eigenstates of the DQD is needed to obtain the few-particle scattering 
states. To this aim, a time-independent approach based on the QTBM has been 
used [2TJ [221 123] . Its solution gives the reflection and transmission amplitude of each 
scattering channel as a function of the initial energy of the incoming electron. Such an 
approach permits to analyze the relation between the resonances in the transmission and 
reflection spectra and the appearance of quantum correlations between the particles, as 
already pointed out elsewhere [25], |23} [26] . All the travelling components of the scattered 
carrier, both reflected and transmitted, have been used to evaluate the creation of the 
entanglement between the dots, together with their decoherence. 

Our numerical simulations show that as a consequence of the scattering between 
an electron injected with a suitable energy and two electrons bound in the ground 
state of the DQD system, the latter can be excited, ending up in an entangled state 
of the constituent dots. This process leads to the appearance of resonance peaks and 
dips in transmission and reflection spectra of the first-excited and ground scattering 
channels, respectively. A side-effect of such a scattering is the loss of quantum coherence 
of the DQD as whole due to its coupling to the scattered carrier. The condition of 
maximum entanglement between the two dots is reached when the bound electrons are 
fully raised to the two-particle first-excited level of the DQD system (which corresponds 
to a Bell-state formed with the single-particle ground and first-excited states of the 
two dots). In this case, the DQD decoherence is zero, since a single output channel 
is possible. However, a single collision is not able to fully excite the dots. We found 
that, in order to build up the maximum amount of quantum correlation between them, 
a repeated injection of charge carriers, that is an electric current, is needed. Indeed, 
at each scattering event the excitation probability of the dots increases until it reaches 
asymptotically 1, which means that a Bell state is obtained, fully decoupled from the 
degrees of freedom of the scattered carriers. In other words, the Coulomb interaction 
between an electric current and two electrons bound in the ground state of a DQD 
structure allows for the maximum entanglement production while the decoherence effects 
on the system vanish. This is in agreement with the procedures adopted in other 
works [261 El SHI ESI EQ], where the entangling schemes are based on the successive 
interactions of a mediator with the qubits. However, in our scheme, the indirect coupling 
of the two dots due to interaction with the scattered carriers, can produce disentangling 
effects, as well. Indeed, a proper tuning of the electric current makes the DQD, initially 
in a Bell state to relax into the ground state, with no quantum correlations. Also in 
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this case the process results to be robust against decoherence. 

Finally, the results here reported show how a suitable electron current, where all 
the carriers have almost a given kinetic energy, permits to switch coherently on and off 
the entanglement between the dots of a DQD structure. Although several interaction 
mechanisms, as electron-phonon coupling can lead to the loss of quantum coherence of 
the DQD in a real experimental setup, we showed that interaction with the mediator 
electrons does not generate entanglement with the leads. Thus, no intrinsic decoherence 
is implied. 

Appendix A. Evaluation of the reduced density matrix 

Here we shall give the explicit derivation of the reduced density matrix of the two dots 
(see equation (TlTl)). initially taken in their ground state, once n carriers injected with 
kinetic energy To close to Tq have been scattered. In order to simplify the calculation, 
the basis C = {|e ), |ei), |e 2 ), |e 3 )} of the DQD eigenstates will be used. Once obtained 
the reduced density operator pr in the C basis, its expression p^ 1 in terms of the E is 
straightforward (see Table [[]) . 

The output three-particle state of equation f|T6|) . stemming from the scattering 
between one carrier injected in the device with Tq=Tq and the two electrons in the 
ground state of the DQD, can be written as 

I^ut) = (boo\T< (1) ) + Qx)|T > (1) >) |e ) + (%, 2 |T< 2 (1) > + 5 02 |T_> (1) >) |e 2 >, (A.l) 

where the superscript (1) means 1 carrier injected, and the reduced density matrix of 
the DQD system can be obtained by tracing y& OU T /y&out I over the degrees of freedom 
of the scattered carrier 
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When a second carrier is injected, after the exit of the previous one from the scattering 
region, the new output state is 

l*&V> = (b o\T< {2) ) + Qx)|T >(2) >) (6oo|T <(1) ) +5 o|T >(1) )) |e > + 

(%2|r< (2) > +c 02 |T> (2) >) (6oo|T< (1) > + qx,|T> (1) >) |e 2 ) + 

(W <(2) > + £ 22 |T> (2) >) (fe 02 |T< 2 (1) ) + £o2|T> (1) >) |e 2 >. (A.3) 

As stressed in Section HI when the DQD is in |e 2 ), the scattering event does not produce 
the relaxation of the dots which remain in the excited energy level. This means that 
in the above expression |6 22 | 2 + |c 22 | 2 = 1. The reduced density matrix of the DQD 
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where Poo=|^oo| 2 + |coo| 2 =l _ | &02 1 2 — | c 2 1 2 have been used. For the case of n scattered 
particles we get 
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(A.5) 



as derived by induction in the following. Assume that expression (IA.5j) is true for n. 
This implies that p^ 1 = Pool e o)( e o| + (1 — Poo) I e s) ( e 2 1 • After the injection of the (n + 1)- 
th carrier, the density matrix p'(ei, e v T^ >y \ ^) describing the total system 

can be evaluated from : 



<(>)(«+!) 



n+1) 



Poo 



b 00 \T< {n+1) ) +£oo|T > ( " +1) )) |6 )(6o| (b* 00 (T<^\ + £ 00 <T > (n+1) 
' + 5oo|T > {n+1) )) |e )(e 2 | (b* 02 (T< + c* 02 <T> 
hr 2 \T< 2 ^) + 5 02 |Tr 2 (n+1 M |e 2 )(e | (W< (n ^\ + c* (T >- ' | ) + 
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,>(n+l) 



~b 02 \T< 2 + c 02 |T> |e 2 )(e 2 | fe(T< 



(1-/4) |^|^ (n+1) ) +£ 22 |^ (n+1 M |6 2 }(e 2 | fe(T< ( " +1) | + c 22 (T > ( " +1) 



> (n+1) 



+ 

(A.6) 

By tracing p'(ei, 6j, \Tm^ >y ) over the degrees of freedom T^ >y ^ of the 

carrier, one obtains the reduced density matrix pi n+1 ^ of the DQD scattered by n + 1 
electrons. Its expression reads 

Pl n+1) ' = Poo(M 2 + Icool 2 ) |6 }(e | + (Ko(M 2 + |c 02 | 2 ) + (1 - p n 00 )(\h 2 \ 2 + |c 22 | 2 )) |e 2 )(e 2 | = 

Poo +1 ko)(e | + (l-Poo +1 )k 2 )(6 2 |. (A.7) 

Thus the expression flA.5j) is true for n + 1. 

Finally, the unitary transformation of Table [1] can be applied to p^ in order to 
obtain the reduced density matrix of the two electrons bound in the DQD given in 
equation ( TT7|) . 
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